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Multiplicative logarithmic corrections frequently characterize critical behaviour in statistical 
physics. Here, a recently proposed theory relating the exponents of such terms is extended to 
account for circumstances which often occur when the leading specific-heat critical exponent van- 
ishes. Also, the theory is widened to encompass the correlation function. The new relations are then 
confronted with results from the literature and some new predictions for logarithmic corrections in 
certain models are made. 
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In a recent Letter, we presented three new relations be- 
tween some of the exponents of multiplicative logarithmic 
corrections to scaling which are frequently manifest in 
statistical physics 0. While two of these relations were 
demonstrated to hold generally, the third fails in certain 
circumstances where the specific-heat leading exponent 
a vanishes. Here, a broad theory which encompasses 
such scenarios is presented and a fourth general scaling 
relation for such logarithms is derived. Together, these 
amount to logarithmic analogs of the standard scaling 
relations between the leading exponents, which are well 
established and fundamentally important in statistical 
mechanics (see e.g. Refs. [2j,|3| and references therein). 

Denoting the reduced temperature by t, we address 
the circumstances in d dimensions where the correlation 
length, specific heat, susceptibility and correlation func- 
tion scale respectively as 
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in which x represents position on the lattice, whose extent 
is indicated by the subscript. When this is finite, the 
counterpart of QJ may be assumed to be Q 



£t(0)~L(lnL)*. 
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The aim of this Letter is to establish full logarithmic 
analogs of the following two standard scaling relations; 



vd = 2 — a , 
v(2-r]) = 7. 



(6) 
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The relation © was developed by Widom jj| (see also 
Ref. 0) who also showed how a logarithmic singularity 
may arise in the specific heat if a = with, in general, a 



superimposed finite discontinuity (see also Ref. [6j). The 
second relation J7J) is due to Fisher |7j . For an authorita- 
tive and comprehensive outline of the development of the 
original scaling relations the reader is referred to Ref. |2j • 

The scaling theory presented herein is based on self- 
consistencies, which are manifest as relations between 
the various correction exponents. For ab inito model- 
specific theories, the renormalization group and related 
approaches are appropriate Q and the reader is again 
referred to Ref. [2j for a review. Our approach is not de- 
pendent on such renormalization group considerations. 

In Ref. 0, we used a Lee- Yang analysis to establish 
the following scaling relation between the exponents of 
the logarithmic corrections analogous to 



a = dq — dv . 



(8) 



and this formula was confronted with a variety of results 
from the literature. While it holds in most models, ex- 
ceptional cases that were identified include the pure Ising 
model in two dimensions and its uncorrelated, quenched, 
random disordered counterpart. Indeed, it is not alto- 
gether surprising that a Lee- Yang analysis, which focuses 
on a complex odd (magnetic) scaling field, cannot com- 
pletely realize the general relationship between the even 
correction exponents appearing in (JHJ). The first main 
aim of this Letter is to redress this situation by appeal- 
ing to Fisher zeros, which are appropriate to the even sec- 
tor and to present a complete theory for the logarithmic 
analog of JBJ), which also neatly encapsulates the d = 2 
(pure and random) Ising cases. The second main aim 
is to present a complementary analog of Fisher's scaling 
relation J2J involving the correlation-function correction 
exponent f] and confront it with the literature. We now 
address these two issues in sequence. 

In two dimensions, the pure Ising model has a = 1 and, 
since all other correction exponents vanish, (JHJ fails there. 
It also fails in the version with random-bond disorder, 
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where q = H03, d = 0, z> = 1/2 and [U03 

~ln|ln|t||. (9) 

Numerical works supportive of the vanishing of d and the 
double-logarithmically divergent specific heat are found 

in Refs. [1 El El ( see also Ref - E3) and Ref - E3 for 
both the random-bond and site-diluted models, respec- 
tively. However, there has been considerable disparity 
in the literature as to the precise scaling behavior and 
counter claims that the specific heat remains finite in the 
random-bond E3> 0] and random-site versions 0, H3| 
also exist (see also Ref. [2ljv 

While it was mooted in Ref. that the detailed loga- 
rithmic corrections in the random-bond and random-site 
Ising models in d = 2 dimensions may in principle differ, 
it is herein clarified that this is not, in fact, expected to 
be the case. Whereas in Ref. |lj , the Lee- Yang zeros were 
used to link the even and odd scaling fields, we now ap- 
peal to the Fisher zeros of the even sector [6j, as that is 
where the apparent specific heat anomaly related to (jHJ) 
lies. The puzzle is resolved as being due to two special 
properties of the pure and random Ising models, namely 
the vanishing of a and the manner in which the Fisher 
zeros in these models impact onto the real axis. 

From the finite-size scaling (FSS) hypothesis, one has, 
for the specific heat |2^ . 



Cl(0) 

Coo(t) 



q(o) 
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Fixing the scaling ratio £l (0)/£oo (t) gives t 
L-V" (In which from @ yields 

C L (0) ~ L^(\nLf- a ^ . 



(11) 



A FSS theory for partition function zeros for pure power- 
law scaling was formulated in Ref. 23] by writing the 
partition function for a finite-size system as a function 
of the scaling ratio there. Here, allowing for logarith- 
mic corrections, this partition function may be written as 
Zh{t) — Q (£l(0)/£oo(*)) and vanishes at a Fisher zero. 
Labeling the j th such zero as tj (L) , one has 



q(Q) 



= Q7 1 (0). 



(12) 



where Q 1 (0) is the j th complex root of Q. Therefore 



\t j (L)\~L--(\nL)- 



(13) 



No assumptions other than the validity of FSS have been 
used to derive (fTT|) and (fTT^) . 

The total number of conjugate pairs of zeros, A/", in a 
suitable variable t is proportional to the lattice volume 
so that Af oc L d . The full expression for the scaling of 
the j th zero is given in Ref. Q (see also Refs. [H E^) 



as a function of a fraction of the total number of zeros 
(2j — 1) /2L d . Then, allowing for logarithmic corrections, 
(f 1 31) is more appropriately written as 



J - 1/2 
L d 



In f :*-_/_) j exp (#;(£)), 

(14) 

where 4>j(L) is the argument of the jth zero. In all 
known cases, the Fisher zeros for isotropic models on ho- 
mopolygonal lattices lie on curves in the complex plane 
and impact onto the real axis along a singular line |26| . 
We assume this scenario, and denote the impact angle 
onto the real axis in the thermodynamic limit by 4>. 

Now, writing the finite-size partition function in terms 
of its Fisher zeros or free-energy singularities, 
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where tj(L) and t*{L) are complex conjugate pairs. As- 
sume that the A4 oc Af zeros which dominate scaling 
behavior close to the critical point are described by the 
scaling form i|14[l . Appropriate differentiation gives for 
the specific heat at t = 



M 



C L (0)~ -L- d ReJ2tJ 2 (L) 



(16) 



.7 = 1 



having included the volume factor L d . 

In the case where vd ^ 2, so that a ^ by ©, one 
finds that the FSS expression (|16fl gives for the singular 
part of the specific heat 



Cl(0) 



(InL) 



(17) 



Comparison with ljll|l leads to the recovery of JHJ . 

If, however, a — 0, so that vd = 2 by ©, the FSS 
expression (llfit for the specific heat becomes 



Cl(0) 



M 

E 



cos(2^(L)) 
J - 1/2 



In 



J - 1/2 
L d 



(18) 



For sufficiently large L and close to the transition point, 
<pj(L) ~ <p and the cosine term in l|18fl becomes a non- 
zero constant provided <f> ^ 7r/4. This is the case in 
the square-lattice pure Ising model in d — 2 dimensions, 
where <f> — ir/2 0. Simple invariance symmetries (such 
as self-duality or duality combined with the star-triangle 
relation) which the distribution of Fisher zeros must re- 
spect ensure that this is also the case with the pure 
model on other lattices [27l| as well as for the symmetric 
random-bond counterpart [2j|. On continuity grounds, 
one also expects 4> ^ 7r/4 in the general random-bond 
and random-site Ising models in two dimensions. 
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In these cases, from the Euler-Maclaurin formula, the 
leading scaling behavior for large L when a — is 



Cl(0) 



(InL) 1 - 21 
lnlnL 



if 2{v -q) + v 
if 2(t> -q)=v 



(19) 



In the thermodyamic limit, \t\ and Coo{t) replace L and 
Cl(0) in l(T9^l. respectively. Comparing JHJ with i fTHjl 
and using ©, one finds 



a = 1 + dq — dv . 



(20) 



This formula replaces (JSJ in such circumstances where the 
model has a — and <f> ^ 7r/4. In the pure Ising model in 
<i = 2 dimensions, where q = v = 0, (|2f)|> gives, correctly, 
q = 1 and the divergence of the specific heat there is 
caused by the extra logarithm as compared with © (see 
also Refs. 0, 0, 01)- I n the random d = 2 Ising model 
where q = 0, v — 1/2, it gives a = 0. In general, if 2{v — 
q) = v, the specific heat instead diverges with a double 
logarithm after l|19|) ■ This is precisely the circumstances 
in the random Ising model in two dimensions [TH . 

The A-colour Ashkin- Teller model also has a = and 
is self-dual El with a = -N/(N - 2) and v = (N - 
1)/(A — 2) 11]. If q = 0, these values also support the 
new scaling relation Q20JI. 

The O(N) symmetric 4> 4 theories (with short- or long- 
range interactions) at their upper critical dimension also 
have a = 0. There, however, (j> = tt/A HHHH^ so that 
(121 H does not follow from H18[l . Instead JHJ remains valid 
there as demonstrated in Ref. Q. 

We now turn our attention to the correlation function 
Q and a new scaling relation for fj, analogous to 10. 
Firstly, fixing the argument of the function D in Q , 
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Following Ref. |3j for example, and writing the singular 
part of the magnetic susceptibility as 



one obtains 



f£oo(*) 

Xoo(*) = / d d xQ 00 (x,t) . 
Jo 
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(22) 



(23) 



From and Q), the leading scaling recovers J7J. Match- 
ing the logarithmic corrections yields 



V = 7 - H 2 - v) 



(24) 



This approach, obtaining the susceptibility from the cor- 
relation function, comes from the original one used by 
Fisher and has also been used in Ref. [13] for the 
d = 2 four-state Potts model. In fact, there rj = 1/4, 
7 = 3/4, v = 1/2, 77 = -1/8 Sill and holds. 



For average quantities in the random Ising models in 
d = 2 dimensions, 77 = 1/4, 7 = 7/8, v = 1/2 and fj = 
G3 and (|21 is again obeyed. This value for_n has been 

convincingly verified numerically mini in mm. The 

new relation l|24l) also holds in the iV-colour Ashkin- Teller 
model, which, along with 77 = 1/4 7 = 7(N-l)/4(N-2), 
v = {N ~ 1)/(N - 2), has 77 = (lj. 

The O(N) symmetric 4 theories at their upper critical 
dimension d = d c = A have 77 = 7 = (AT + 2)/(A r + 8), 
v = (N + 2)/2(N + 8) and 77 = 34], and the expression 
(I24|l is obeyed. Likewise, O(N) spin models with long- 
range interactions decaying as x _ ^ +<J ' ) have logarithmic 
corrections at d — d c — 2a. There, 77 = 2 — a [35|. 
7 = (A + 2)/(A + 8), v = (N + 2Vo-(A + 8) and the 
relation (|24l) correctly yields 77 = |3q . 

For the percolation problem, 77 = 0, 7 = 2/7 , and 
v = 5/42 at the upper critical dimension d c = 6 |37| . The 
correction exponent for the correlation function there has 
recently been calculated to be 77 = 1/21 j3^. Again, this 
set of values satisfies (|24|) . 

Finally, (|24|l can be used to predict the value of 77 
in other models, such as m-component spin glasses and 
Yang-Lee edge problems at their upper critical dimension 
d c = 6. For the former, 77 = |39j, 7 = 2m/(2m — 1), 
v = 5ra/6(2m— 1) 0] giving 77 = m/3(2m— lL For the 
Yang-Lee problem, 77 = 0, 7 — 2/3, v = 5/18 |40j so that 
the prediction from H24[) is 77 = 1/9. These values remain 
to be verified numerically. 

It is observed in Ref. that the magnetization 

77i 00 (/j) for the four-state Potts model may be deduced 
from the correlation function by an alternative argument; 
representing a generic spin- type variable by s(x), if the 
spins decorrelate in the limit where x — > 00 such that 
SooOM) = (s(O)s(z)) -» (s(0))(s») = m^(t) there, 
then using and (|21|l and matching with m 00 (<) ~ 
|*|^| In |t||^ Q gives i/(d-2 + r?) = 2/3 and 7} = 2/3 + £>(d- 
2 + 77). From the standard scaling relations, the first of 
these again recovers J7J). (See also Ref. |41|.) From the 
scaling relations for logarithmic correction [l| the second 
yields 7} = dq + 7 — z>(2 — 77). When g vanishes, this is 
identical to (|24|) . Indeed, this is the case in the d = 2, 
four-state Potts model [lias well as in the d = 2 pure 
and random Ising modelsfll Wi l32| . However, since q ^ 
at the upper critical dimension of the O(N) <fi 4 theories 
and their long-range counterparts, the percolation prob- 
lem, spin glasses and the Yang-Lee problem 0, |4(J, |42| , 
this detailed matching of ^00(2;, t) with m^ Q (/j) is invalid 
in these cases. Instead, l|2*4^) holds in each case. 

In conclusion, then, the scaling theory presented in 
Ref. m nas been extended to deal with the specific heat 
when its leading exponent a vanishes and the Fisher ze- 
ros impact onto the real axis at an angle other than it /A. 
In such cases, (JHJ is replaced by pHjl. Also, the gen- 
eral theory has been extended to deal with the correla- 
tion function and the new relation (|24|l has been checked 
against the literature and predictions made. Together 
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with Ref. , the new formulae l(2Tj)l and l(24*|l offer a set 
of scaling relations analogous to the standard ones and 
appropriate to logarithmic corrections. 

Besides these general results, progress specific to the 
random Ising models in two dimensions has been made. 
Through l|2UI) and (|24[) . the hitherto numerically most 
elusive and controversial quantity a has been directly re- 
lated to rj, fj, 7 and q, all of which are clearly established. 
Moreover, our theory automatically generates the famous 
double logarithm in the specific heat in these instances. 
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